Deep Strong Coupling Regime of the Jaynes-Cummings model 
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We study the quantum dynamics of a two-level system interacting with a quantized harmonic oscillator in the 
deep strong coupling regime (DSC) of the Jaynes-Cummings model, that is, when the coupling strength g is 
comparable or larger than the oscillator frequency u (g/oo > 1). In this case, the rotating- wave approximation 
cannot be applied or treated perturbatively in general. We propose an intuitive and predictive physical frame to 
describe the DSC regime where photon number wavepackets bounce back and forth along parity chains of the 
Hilbert space, while producing collapse and revivals of the initial population. We exemplify our physical frame 
with numerical and analytical considerations in the qubit population, photon statistics, and Wigner phase space. 



The interaction between a two-level system and a harmonic 
oscillator is ubiquitous in different physical setups, ranging 
from quantum optics to condensed matter and applications to 
quantum information. Typically, due to the parameter acces- 
sibility of most experiments, the rotating- wave approximation 
(RWA) can be applied producing a solvable dynamics called 
the Jaynes-Cummings (JC) model Q. In this case, Rabi os- 
cillations inside the JC doublets or collapses and revivals of 
the system populations |2] are paradigmatic examples of the 
intuitive physics behind the JC dynamics. To achieve these 
and other phenomena in the lab, the strong coupling (SC) 
regime is required, that is, the qubit-oscillator coupling has 
to be comparable or larger than all decoherence rates. This 
model accurately describes the dynamics of cavity QED 0[4]|, 
trapped ion experiments [5 ], and several setups in mesoscopic 
physics, where the qubit-oscillator model is essential in mod- 
eling superconducting qubits (6) with either coplanar trans- 
mission lines (TUTOl or nanomechanical resonators lfTTl[T2ll . 
Nowadays, solid-state semiconductor [13 ] or superconductor 
systems |[T4lj20l have allowed the advent of the ultrastrong 
coupling (USC) regime, where the coupling strength is com- 
parable or larger than appreciable fractions of the mode fre- 
quency: g/uj > 0.1. In this regime, the RWA breaks down and 
the model becomes analytically unsolvable, although some 
limits can be explored f[2TI - [25l . Confident of the impressive 
fast development of current technology, one could explore fur- 
ther regimes where the rate between the coupling strength and 
oscillator frequency could reach g/uj > 1, here called deep 
SC (DSC) regime. This unusual regime, yet to be experimen- 
tally explored, is the focus of our current efforts. In this letter, 
we introduce a rigorous and intuitive description of the DSC 
regime of the JC model, providing an insightful picture where 
photon number wavepackets propagate coherently along two 
independent parity chains of states. In this way, the Hilbert 
space splits in two independent chains, exhibiting a compre- 
hensible collapse-revival pattern of the system populations. 

We consider the Jaynes-Cummings Hamiltonian without 
the RWA, also called Rabi Hamiltonian, describing a two- 
level system coupled to a single mode harmonic oscillator 

H = —(jJocr z + Tiuja) a + hg(a + + a~)(a + a^). (1) 



Here, a and a) are the annihilation and creation operators of 
the mode with frequency cj, while a z and a ± are Pauli oper- 
ators associated to a qubit with ground state |g), excited state 
|e), and transition frequency uq. We concentrate in the study 
of the DSC regime, g/uo > 1, with no particular relation be- 
tween uj and ujo. We do not refer to any particular system be- 
cause several of them in quantum optics and condensed matter 
may profit from the physical insight developed here l26l . We 
start by observing that the parity operator l27l 

n = -<x z (-ir° = -(|e)(e| - |g)(g|)(-l) ot °, (2) 

with U\p) = p\p) and p = ±1, is a key element associated 
to the Hamiltonian in Eq. ([I]). It is instrumental to understand 
how the system dynamics moves inside the Hilbert space split 
in two unconnected subspaces or parity chains, 

|g0 a > |el a > |g2 a > |e3 a ) . . . (p = +1), 
|e0 a ) <-> |gl a > O |e2 a > |g3 a ) . . . (p = -1). (3) 

Neighboring states within each parity chain may be connected 
via either rotating or counter-rotating terms. For example, 
in the parity chain with p = +1, the counter-rotating term 
cr + a^ induces the transition |g2 a ) — > |e3 a ), while the rotating 
term cr + a induces |el a ) |g2 a ). When going back from 
DSC^USC^SC, the parity chains break into the known 
Jaynes-Cummings doublets {|g, n a + 1), |e, n a )} because we 
enter into the domain of applicability of the RWA. 

We introduce the parity basis \p,rib), where b^b\rib) = 
nb\nb), and b = cr x a such that b\p,rib) = y/nb\p, rib — 1). 
Using this basis, the Hamiltonian in Eq. ([T]) can be written as 

H = hutfb + hg(b + 6 f ) - fry (-l) 6t6 n. (4) 

This Hamiltonian commutes with the parity operator n, and 
for each parity chain (p = ±1) there is an independent Hamil- 
tonian describing a perturbed harmonic oscillator. Note that 
the term — Tiujq(— l) 6t6 n/2 behaves as an energy shift pro- 
portional to uq. In the DSC regime, we can get rid of the term 
hg(b + fot) in Eq. ^ by changing to the basis D(—f3o) |p, rib), 
with D(f3o) = e^° b and /3q = g/uj. The eigenenergies 
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FIG. 1. (Color online) (a)-(b) Round trip of a photon number 
wavepacket and collapse-revivals due to DSC dynamics with initial 
state | + ,0&) = |g, a ). c) Collapse-revivals with secondary peaks 
due to counterpropagating photon number wavepackets starting in 
initial state |+, 2b) — |g, 2 a ). For all cases, cjq = and g/uo — 2. 



and eigenfunctions can be approximated as a series in uq/uj 



Efajh « ujn h - g 2 /oj - ^ p (-l)^A nbnb + 
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Alternative approximations can be found in the literature l2T1l . 
To first order we get a displacement in the energy levels due 
to the coupling A nbUb = (nb\D(2/3o)\nb), a correction which 
is much smaller than one, |A nbm J <C 2~ (nb+mb) . Note that 
this formalism is rigorously valid in the DSC regime. 

We study now the DSC dynamics with the initial state 
\ip(0)) = | +,05) = |g, a ), as we activate the interaction 
in Eq. We observe that the photon statistics, P nb (t), 
will spread independently along each parity chain, eventually 
reaching an energy barrier and bouncing repeatedly. Remark- 
ably, an intuitive picture can be found, as displayed in Figs.[T] 
and [2] , that provides physical insight in a problem that is, in 
general, analytically intractable. Note that, in Figs. [T^ and 
[TJ), the round trip of the initial photon number wave packet 
induces collapse-revivals that are not reminiscent of the SC 
regime of the JC model [ 2 ], where initial large coherent states 
are required. In the DSC limit, with uo = 0, this intuitive 
picture can be rigorously confirmed integrating the evolution 



|V(*)> = D t (/?o)e- i(wbtb - s2/w) *D(/? )|+,O fc ) 
= tf(t ) w = 0)|V(0)> = e 



(6) 



where f3(t) = (3o(e ZUJt — 1) is the amplitude of a coherent 
state. The revival probability of the initial state reads 



p +0b (t) = \(m\m)f 
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FIG. 2. (Color online) (a) Photon statistics at different times of the 
evolution with c^o = 0.5 uj. (b) Comparison of probability P+,o b if) 
calculated for cjo = (solid line) and cjo = 0.5 cj (dashed line). In 
all simulations the initial state is |+, 0&) and g/u — 2. 



exhibiting periodic collapses and full revivals (28). When 
the initial state is |+,2&) = |g, 2 a ), as in Fig. the 
DSC dynamics generates counterpropagating photon number 
wavepackets in both directions that bounce back and forth 
producing interference secondary peaks. Similar intuition 
follows when considering initial superposition states, e. g. 
(|+, O5) + |+, 2fr))/\/2, as long as the state components be- 
long to the same parity chain, otherwise no secondary peaks 
appear. When we break the qubit degeneracy, ujq ^ 0, the 
intuitive picture remains but we lose the integrability of the 
problem. Probability still spreads along each parity chain, as 
seen in Fig. [2] but now the photon number wavepacket suffers 
self-interference, it distorts and its center no longer follows 
the periodic orbits of ujq = 0. The result are full collapses and 
partial revivals where probability P + o b is not completely re- 
stored, and whose maximum value deteriorate as time passes. 

The collapses and revivals have also interesting conse- 
quences in phase space, which we have analyzed using the 
Wigner function and phase space trajectories, (x(t),p(t)) = 
((b + ib ] - ib)/V2). In the integrable case, ujq = 0, the 
Wigner function of the state is a Gaussian centered on a point 
x + ip = j3(t) which draws periodic circular orbits on the 
plane. As soon as we switch on the term proportional to c^o, 
the wavepacket suffers two distortions, see Fig. [3] The first 
one is a squeezing tangential to the orbit, shown in Fig. |3j:. 
Accompanied by the difussion and interference in the Wigner 
wavepacket, the orbits also distort, becoming spirals that relax 
towards the center of the original orbits, — /3q. 

The phenomenon of collapses and revivals for the noninte- 
grable case, ujq ^ 0, even if partial, reveal a structure in the 
Hamiltonian spectrum, which is approximately equispaced. 
We can write the revival probability 



F +06 (i) = |^|^(0)|^)| 2 e-^/ ft | 2 , 



(8) 



(7) as a function of the overlap of the initial state with the eigen- 
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FIG. 3. (Color online) Modulus of the Wigner function, |W(x,p)| 5 
and trajectory (x, _p), for Eq. (a) For ujo = the Wigner function 
is a symmetric Gaussian (a coherent state) and moves clockwise in 
phase space along a circle. (b)-(c)-(d) For ojo = 0.5 cj the Wigner 
function loses its symmetry and moves clockwise along a spiral tra- 
jectory as shown in the sequence corresponding to t — 0.5, 1,5. In 
both simulations the initial state is | + , 0*,) and g/oj = 2. 
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FIG. 4. (Color online) Collapses and revivals of P+,o b (t) for g/oj — 
2 and (a) ujq — 0.3cj and (b) ujq — 0.5cj. We plot the exact numerical 
solution (area) at first order in ujq /oj (solid) and a two-mode approx- 
imation (red, dashed), (c) Distribution of probability of the different 
detunings 5 nb = (Tiujrib — E nb ) / Tiujq , weighted by their contribution 
to the wavefuncion given in Eq. (8]). We have marked the level N£ 
which is used for the two-mode approximation in curves (a) and (b). 



states of the full model, (H - Ei)\<j>i) = 0. When uj = 0, the 
eigenenergies are regularly spaced, Eg — huj£ and the func- 
tion becomes periodic with period 2tt/oj. This causes an ini- 
tial Gaussian wavepacket in phase space to get reconstructed 

at the same position for t = 27r, 47r, 67r, In the DSC case 

with ujq ^ 0, the energy levels deviate very slightly from this 
regular distribution, Eg = huj£ — hSgoJo, where the correc- 
tion 5gujQ is less than 10% for the examples considered in 
this work. The reconstruction of the wavepackets is incom- 
plete and different partial waves may get delayed or acceler- 
ated with respect to the original orbit, f3(t). This causes the 
squeezing of the Wigner function and the self-interference in 
the photon number wavepacket, as displayed in Fig. [2^. 

We have also found that the overall dynamics is very ac- 
curately captured by the first order correction to the eigenen- 
ergies, shown in Eq. If we use 5 Ub = ( — l) nb pA nbTlb /2 
in Eq. together with the initial condition (ip(0)\(j) nb ) = 
(1>(0)\D(-Po)\pn b ) = 5^ p exp(-\(3\ 2 /2)(3^/n b l we ob- 
tain curves that approximate very well the exact result. This 
is shown in Figs. [4^ and[4|) for values of ujq = 0.3uj and 0.5a; 
in the DSC regime, g/uj = 2. In both cases the revivals hap- 
pen close to t = k (2ir/u), with integer k, but decreasing in 
intensity and with a large fraction of the curve moving with a 
slower speed, reconstructing itself at later times. The effect of 
this is more evident in Fig. [4]}, and also in the photon number 
wavepacket plot in Fig. [3^, where one appreciates two waves 
with slightly different periods interfering with each other. The 
existence of this delayed revivals is due to the structure of the 
wavefunction ip(t) that, as shown in Fig. [3J:, is composed of 
many contributions close to zero detuning, S nb ~ 0, and a 
few large contributions with S = —0.116,-0.223. The for- 
mer constitute the main revivals, while the second ones make 



the revivals at slightly longer periods, 2ttuj— 0.223cjo, forming 
the second wavefront in Fig. [2^ and Fig.[4^-b. 

Based on our previous results, we have developed a heuris- 
tic approximation that allows us to reproduce the main re- 
vivals. Our method recognizes that if we start the dynam- 
ics with state ^(O)) = |+, iV&), the main contribution to the 
wavefunction is around a level N£ = [(g/uj) 2 + A/5], where 
[•] denotes the closest integer. This is indeed the case in the 
considered examples, as shown in Fig. [4]:. We will only con- 
sider the energy correction for this level, S^r , and neglect the 
dephasing of all other off- re sonant terms. Under this criteria, 
we approximate the system state, up to normalization, as 

\m) «i/(t,^o = o)i^(o)> +^r( C - <(wJVb - woA ^^ /2)t 



-^)D(-(3 )\+,N r b ), 



(9) 



where U(t, ojq = 0) is the evolution operator at ujo = intro- 
duced in Eq. ([6]), and ifjjyr = (<p N r\^(0)). This form of the 
state is motivated by the behavior of the Wigner function, as 
seen in Fig.[3p3): a central core approximated by the solution 
of ujo — plus a delayed correction capturing the effects of 
co>o, forming the squeezed tail. Considering the simplest case 
of Nb = 0, from state in Eq. ([9]), it is straightforward to obtain 
a simple analytical expression for the revival probability 



-\{3(t)\ 2 /2-(3* 



-IW)I 2 



,2iVT 



cos(co>o^/v r £/2) — 1 



(10) 



where small terms are neglected. This expression has been 
compared with the exact solution, as shown in Fig. [4^ and 
[4|), giving a good estimate of the height of the partial revivals 
as a function of time and ujq. The other features, such as the 
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FIG. 5. (Color online) Parity chains for the two-qubit DSC regime. 



delayed front is not reproduced because this approximation 
does not contain the contributions with S = —0.223, but this 
can be improved by including more resonant levels. 

Two qubits and a mode.— It is also possible to give a quali- 
tative and quantitative description of the DSC considering the 
case of two qubits. Here, the 2 (g) 2 (g) TV-dimensional Hilbert 
space bifurcates into two independent parity chains of tetra- 
hedra, see Fig. |5j where each vertex is connected to their 
neighbors via rotating or counter-rotating terms. The same 
dynamical properties of probability collapses and revivals can 
be found, as well as interesting entanglement properties. 

Conclusions.— The SC regime of the JC model is consid- 
ered nowadays an intuitive and comprehensive field. The USC 
regime is described by the SC regime plus RWA and higher- 
order corrections. In this work, we have aimed at developing 
an insightful description of the DSC regime of the JC model. 
The transition between the SC and DSC regimes remains a 
rather diffuse crossover with well understood frontiers. 
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